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ABSTRAGT 

We have observed and spatially resolved a set of seven A-type stars in the nearby Ursa Major moving 
group with the Glassic, GLIMB, and PAVO beam combiners on the CHARA Array. At least four of 
these stars have large rotational velocities (iisinz > 170 km s“®) and are expected to be oblate. These 
interferometric measurements, the stars’ observed photometric energy distributions, and wsinz values 
are used to computationally construct model oblate stars from which stellar properties (inclination, 
rotational velocity, and the radius and effective temperature as a functi on of latitude, etc.) are 
determined. The results are compared with MESA stellar evolution models (|Paxton et al.|2011 2013| 
to determine masses and ages. The value of this new technique is that it enables the estiniation of the 
fundamental properties of rapidly rotating stars without the need to fully image the star. It can thus 
be applied to stars with sizes comparable to the interferometric resolution limit as opposed to those 
that are several times larger than the limit. Under the assumption of coevality, the spread in ages 
can be used as a test of both the prescription presented here and the MESA evolutionary code for 
rapidly rotating stars. With our validated technique, we combine these age estimates and determine 
the age of the moving group to be 414 ± 23 Myr, which is consistent with, but much more precise 
than previous estimates. 


Torres et al. 
etimes of 3.6 


1. INTRODUCTION AND MOTIVATION 

Population I stars with spectral type A have masses 
that range from 1.5 - 2.5 Mq, based on dy namical mea¬ 
surem ents of spectroscopic binaries (e.g., 

20101, and corresponding main sequ ence li: 

- 1.0 Gyr (assuming tms M“^'^, Kippenhahn et al. 
20121. However, measuring the age and mass of any 
single A-type star is often severely complicated by their 
peculiar characteristics, including chemically anomalous 
photospheres (e.g., Am stars, Ap stars, A Boo stars), ra¬ 
dial and non-radial pulsations (e.g., 7 Doradus-type, 6 
Scuti-type), and seve rely distorted shapes f rom ra pid ro¬ 
tation (e.g., Reg ulus - McAlister et al.|20'b5 Alta ir - Mon- 
nier et aL]|2007 Vega - |Autdenberg et al.||2006[ |Monnier 
et al. 2012p. Despite these challenges, interest in deter- 
minmg precise ages and masses for A-type stars has been 
heightened with the discoveries of directly imaged candi- 
date planets orbitin g several of t hem, including HR 8799 
(Marois et al.||2008); Fomalhaut (Kalas et al.||2008); etc. 
In all cases the estimated masses of the planetary com¬ 
panions de pend critically on t he ages assigned to the hos t 
stars (e.g., Moya et al. 2010 versus Marois et al. 2010), 
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which in most cases are poorly determined. 

Although a variety of independent techniques exist for 
estimating the ages of Sun-like stars, such as photo- 
spheric lithium depletion or chr omospheric activity (e.g.. 


Mamajek & Hillenbrand 2008), these are ineffective for 
A-type stars with predoinmantly radiative atmospheres. 
For these stars, a more effective method is to compare ob¬ 
servable stellar properties (e.g., radius, luminosity, and 
temperature ) to the predictions of stellar evolutionary 


mode ls (e.g., David & Hillenbrand 2015 Brandt & Huang 


2015). This “isochronal htting” technique has the po- 


tential to work well for A-type stars since their radii, 
temperatures, and luminosities evolve much more sub¬ 
stantially than Sun-like stars do during the first ~Gyr 


of their main sequence lifetime. For example, the MESA 

evolutionary models ( 

Paxton et al. 

2011 

20131) predict 

that the radius, luminosity, and temperai 

ure of a 2 IVIq 


star change by -1-32%, 4-20%, and —10%, respectively, in 
just 500 Myr after the zero-age main sequence}^ This 
can be compared to a 1 Mq star that, in the same time 
frame sees its radius, luminosity, and temperature change 
W only -1-3.0%, 4-8.8%, and 4-0.6%, respectively. Figure 
[b illustrates these differences in evolutionary rates with 
four stars with masses between 1 and 2.5 Mq. 

Unfortunately, the peculiar properties of A-stars make 
these relatively straight-forward comparisons with mod¬ 
els difficult in practice. Pulsation and rapid rotation re¬ 
sult in observable stellar properties that are both time 
and orientation dependent. While the photometric vari¬ 
ations due to pulsations are typically less than a few 


^ The zero-age main sequence is defined for each star to be the 
point at which the contribution to the luminosity of the star due 
to gravitational contraction is ~1% that of core fusion as predicted 
by the MESA evolution code. 
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, rapid rotation distends the 


star so that its size is no longer defined by a single ra¬ 
dius. The resulting gravity darkening that oc curs creates 
a tempe rature gradient on the star’s surface (von Zeipel 
1924ajb|), causing the star to no longer be defined by a 


single temperature. The net effect is that the observed 
flux depends on the star’s inclination, making the total 
luminosity challenging to determine since inclination is 
unknown for most stars. In addition to the challenges 
in observationally determining an A-star’s stellar prop¬ 
erties, their peculiar characteristics must be accounted 
for in the adopted stellar evolution models. Rapid rota¬ 
tion has bee n shown to dramatically affect the way the 
star e volves (Meynet & Maeder|2000 Maeder & Meynet 
20101. For example, the MFSA evolutionary code pre¬ 
dicts that solar- to intermediate-mass stars rotating at 
50% of break-up velocity have average surface tempera¬ 
tures that are significantly cooler relative to non-rotating 
StarBand evolve more slowly off the main sequence (Fig¬ 
ure [H. Finally, the anomalous surface abundances of 
many A-stars can complicate the choice of evolutionary 
model metallicity, which is usually scaled relative to so¬ 
lar. 

Fortunately, with the high angular resolution that opti¬ 
cal/infrared interferometers provide, it is now possible to 


ity darkening laws and stellar evolutionary models that 
include rotation. 

2. THE SAMPLE AND CHARA OBSERVATIONS 

With a nucleus distance of 25 pc, the Ursa Major mov¬ 
ing group is one of the closest and best-studied mov¬ 
ing groups. It consists of 15 nucleus stars and 47 likely 
stream members with an estimat ed age of 500 ± 100 
Myr and a metallicity of Z=0.016 ( King et al.||2003 ). As 
summarized in Table previous studies have touiid an 
age for the moving group ranging from 200 to 1000 Myr . 
The introduction of | Ammler-von Eiff & Guenther (20091 
provides an excellent history of the study of the UMa 
moving group. 

We define a sample of A-stars in the Ursa Major mov¬ 
ing group for interferometric observations by selecting all 
stars with B — V colors less than 0.31 from the “UMa 


nucleus stars” list in King et al. (2003). The hottest of 


these stars, has a B — V color of —0.022 (v an Leeuwen 
2007 ) and an assigned spectral type of Al ( [Gray et al. 


20031. The resulting list consists of 7 stars of which 


2 stars (Mizar A = HD 116656 and Mizar B = HD 
116657) form a spectroscopic binary pair of comparable 
brightness {AMy = 1.68 mag). Mizar A and B are con¬ 
sequently excluded from this sample because the close 


synthesis ( 

Baron et al. 2010 

), to directly determine fun- 

pair would bias interferometric obseryations, making it 

damental 

Droperties of rapidly rotating early-type stars 

difficult to distinguish the physical properties of each star 

(van Belle 

2012 

Che et al. 

2011 

Monnier et al. 

2012 

1 . 

indiyidually. Another of these seyen nucleus stars (Alioth 


In these cases, the oblateness and gravity darkening can 
be observed directly, which enables more accurate deter¬ 
mination of the star’s luminosity and comparisons with 
evolutionary models. However, there are only a handful 
of stars that are large enough and bright enough for this 
technique to work effectively with current facilities. 

We present a technique that allows for the correction of 
the effects of rotational distortion without having to fully 
image the star. Fundamental parameters are determined 
by tuning a model of an observed rapidly rotating star 
such that the model-calculated interferometric visibili¬ 
ties match the observed visibilities obtained at multiple 
baseline orientations; the model is further constrained 
by the star’s photometric energy distribution and pro¬ 
jected rotational velocity (wsinz). The advantage of this 
technique is that it enables the determination of funda¬ 
mental properties of rapidly rotating stars that are too 
small and/or too faint to be observed with imaging in¬ 
terferometric beam combiners. Of the 112 A-type stars 
within 50 parsecs that are observable with the CHARA 
Array (i.e., with S > —10° and with no known compan¬ 
ions within 2" and with AMy < 5 mag), only 13 have 
estimated angular diameters large enough (0^1 mil- 
liarcsecond) to fully benefit from imaging. Another ad¬ 
vantage of this technique is that it does not require the 
measurement of closure phases, so it is not necessary to 
use the many simultaneous baselines that are necessary 
for the aperture synthesis imaging technique. In this 
paper, we demonstrate the success of this technique by 
comparing the relative ages of rapidly and non-rapidly 
rotating stars in the Ursa Major moving group. These 
ages are determined by comparing modeled luminosities 
and radii with the predictions of the MESA evolution 
model. Indirectly, the results thus also provide a new 
age estimate for this moving group and tests of grav- 


(2011) identify a companion to Alioth with a projected 
separation of 0.11" and a AM/ of 2.31 mag. A fourth of 
these seven stars (Alcor = HD 116842) has an observed 
stellar companion of spectral type M3-M4 and with a 


projected separation of 1.11" (Zimmerman et al. 


2010 


Mamajek et al.|2010). However, with a AM// ot ~6, the 


companion is too taint to contaminate the interferomet¬ 
ric observations, so it is not excluded from the sample. 
None of the other nuc leus stars have known companions 
(De Rosa et al. 2014|). The four nucleus member stars 
that are included in this sample are Merak = HD 95418, 
Phecda = HD 103287, Megrez = HD 106591, and Alcor 
= HD 116842. 

There are 6 additional A-stars that are likely stream 
m embers of the mov ing group (listed as “Y” or “Y?” 


King et al. (2003)). Two of these 6 (Menkalinan = 


HD 40183 and Alphecca = HD 139006) are spectroscopic 
binaries with AMy values of ^1 and ~ 4, respectively 


(|Pourbaix||2000| |Tomkin & Popper||1986D and so are not 
observed. Ut the remaining four, one star (21 LMi = 
HD 87696) was not observed due to limited telescope 
time. The remaining three (Chow = HD 141003, 16 Lyr 
= HD 177196, and 59 Dra = HD 180777) are included 
in the sample. One of these stream stars (59 Dra) has a 
candidate brown dwarf companion ( [Galland et al.|20'06 ), 
but this is too faint to contaminate the interferometric 
observations. 

In total, we obtained new interferometric observations 
for 6 Ursa Major A-type stars (3 nuclear members and 
3 stream members). One additional star, Mera k, was 


observed interfe rometrically by a previous study (Boya- 
jian et al. 2012[). These seven stars have spectral types 


rangi ng m)m A0-A7. M erak also has a peculiar metal¬ 
licity (Royer et al.|2014 ) and is an apparent slow rotator 
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(a) (b) 


Figure 1. Plot of temperature versus radius (left) and temperature versus luminosity (right) of the evolution tracks of eight stars with 
masses ranging from 1.0 to 2.5 Mq and an angular rotation rate of either 0 % or 50% that of the break-up velocity. The red circles represent 
the properties of each star while on the zero age main sequence (at 41, 22, 9.5, and 5.7 Myr for the 1.0, 1.5, 2.0 and 2.5 Mq stars, 
respectively for lj = 0.0 and 49, 26, 11, and 5.9 Myr for lj = 0.5) and 500 Myr after that point. 


with a iisinj of 46 ± 2.3 km s“^. While it is possible 
that Merak is a rapidly rotating star oriented pole-on, 
there is some suggestion that the peculiar me tallicity of 
Ap stars is due in part to their slow rotation ( Abt||20d^ 


and references therein). Another apparent slow rotator 
in the observed sample is 59 Dra with a v sin i of 70 ± 3.5 
km s“^. 59 Dra shows a normal A-star metallicity sug¬ 
gesting that it may be a rapidly rotating star oriented 
pole-on. The four stars in this set that are nuclear mem¬ 
bers have distances within the very narrow range of 24.4 
to 25.5 pc, while the three stream members are more 
spread out, having distances of 27.3, 37.4, and 47.6 pc. 
The properties of all seven stars in the set are summa¬ 
rized in Table which includes spectral type, projected 
rotational velocity, Hipparcos distan ce, photometry, an d 
UMa membership as determined by King et al. (2003). 

All observations were obtained using Georgia State 
University’s Center for High Angular Resolution Astron¬ 
omy (CHARA) Array. The CHARA Array is a six 
telescope interferometer w hich operates at optical an d 
near-infrared wavelengths (ten Brummelaar et al.||2005). 
The CHARA Array’s six telescopes are arranged in a 
Y-shaped configuration with baselines ranging from 34- 
331 m. The naming convention for these six telescopes 
consists of a letter representing one of three arms of the 
“Y” (“S” for south, “E” for east, and “W” for west), 
and a number indicating the outer telescope (1) or the 
inner telescope (2) of each arm. Data were obtained us¬ 
ing three beam combiners: Classic, CLIMB, and PAVO. 
All three beam combiners measure the contrast of the 
interference pattern produced by the light from each of 
the telescopes used. This contrast is known as a visi¬ 
bility. The two-telescope Classic beam combiner takes a 
single visibility measurement per observation in a broad¬ 
band near-infrared filter (AT-band for this work). The 
three-telescope CLIMB beam combiner, which also op¬ 
erates in the near-infrared (in either the H- or AT-band), 
takes three simultaneous visibility measurements for each 
broadband observation (one for each combination of two 
telescopes). The PAVO beam combiner was used in its 


two-telescope mode and each observation yields 23 visi¬ 
bilities spectrally dispersed across a wavelengths ranging 
from 0.65-0.79 ^m. Because PAVO and Classic obser¬ 
vations were taken using two telescopes at a time, only 
a narrow range of baseline orientations was used. This 
is illustrated in Figures [2][^ We note that for two stars 
(16 Lyr and 59 Dra), we do not have sufficient baseline 
orientations to measure oblateness. A general observing 
strategy was adopted whereby calibrator stars (described 
in Section were observed both before and after each 
target star. This set of observations is referred to as a 
visibility bracket. Over 8 nights of observing, a total 
of 56 visibility brackets yielding 724 individu al visibil¬ 
ity measurem ents were obtained on 6 stars. [Boyajian 


et al. (|2012 ) obtained 25 brackets on Merak with the 


two-telescope Classic beam combiner. Table lists the 
calibrators, beam combiners, baselines, and wavelengths 
used during each observation as well as how many brack¬ 
ets were obtained for each star. 

3. DATA REDUCTION AND CALIBRATED VISIBILITIES 

Interferometric data from the Classic and CLIMB 
beam combiners were reduced using the reddassic and 


redcli mb pipelines, respectively (Ten Brummelaar et al. 
2013), yielding reduced visibilities tor each observation 


made. The pipeline used to reduce the observatio ns made 
with the PA VO beam combiner is described by Ireland] 
et al. (2008). Many factors, both atmospheric and in¬ 
strumental, serve to decrease the visibility measured by 
an interferometer. This decrease depends in part on at¬ 
mospheric turbulence at the time of observ ation and the 
airmass at which the star is observed (e.g., Boden|[2007 


Roddier 1981). Correcting for these temporal effects on 


the visibility requires frequent observation of a star with 
a known angular diameter that is ideally smaller than the 
interferometric resolution {X/2B). Such a star is called 
a calibrator star. When observed near the target star 
both in time (< 30 minutes) and on the sky (< 10°), 
the target star’s intrinsic visibility (V*) should be ob¬ 
served (VJjj) to be reduced by the same amount as the 
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Table 1 

Age Estimates for the Ursa Major Moving Group. 


Age 
(Myr) 
~300 
SOOitlOO 
630-1000 
300-400 
~500 
~200 
SOOitlOO 
~600 
39 ^^ 
530 ± 40 
414 ± 23 


Reference 




von Hoerner i ^ 

9571 



Giannuzzi i IS 

791 


Eggen 1 199i 

U 


Soderblom et al. 

1993 

1 



Asiain et al. 1 ] 

999 1 



Konig et al. i 2 

0021 

King et al. 1 2( 

)03l 

K 

ing Schuler 

2005 

Id 

a,vid & Hillenbrar 

id (20 

15 1 


Brandt Huang 

(2015 

1 


This work 


^David & Hillenbrand| (|2015|l do not report an age for the UMa^jnoving group. The value listed here corresponds to the median of the 
ages they report tor the 7 Ursa Majoris stars studied here (Table [^. 


Table 2 

Presented Sample. 


Common 

Name 

HD 

Number 

HIP 

Number 

Spectral 

Type“ 

usim^ 

(km/s) 

D<^ 

(pc) 

v4 

(mag) 

B-V^ 

(mag) 

(mag) 

UMa 

Membership^ 

Merak 

95418 

53910 

Al IVps (SrII) 

46 ± 2.3 

24.4 ± 0.1 

2.35 

0.033 

2.285 

Nuclear 

Phecda 

103287 

58001 

Al IV(n) 

178 ± 8.9 

25.5 ± 0.3 

2.43 

0.044 

2.429 

Nuclear 

Megrez 

106591 

59774 

A2 Vn 

233 ± 11.7 

24.7 ± 0.1 

3.34 

0.077 

3.104 

Nuclear 

Alcor 

116842 

65477 

A6 Vnn 

228 ± 11.4 

25.1 ± 0.1 

4.05 

0.169 

3.145 

Nuclear 

Chow 

141003 

77233 

A2 V 

207 ± 10.4 

47.6 ± 0.6 

3.68 

0.073 

3.546 

Stream 

16 Lyr 

177196 

93408 

A7: V 

124 ± 6.2 

37.4 ± 0.2 

5.07 

0.186 

4.505 

Stream 

59 Dra 

180777 

94083 

FO Vs 

70 ± 3.5 t 

27.3 ± 0.1 

5.19 

0.308 

4.313 

Stream 


Notes - (a) Nucleus Stars - G^^et^ ( 2003|, Strea m Stars - ^jevato_fcV^|19|^; (b) ([20071; (c) van Leeuwen 


(20071; (d) [Perryman fc ESA ( |1997| ; (e) (Jutri et al. (2003); (fj [Glebocki fc (jnacinski (20051; (gj King et al. ( 2o63| ). 


calibrator’s (intrinsic - V?, measured - 


W 


V* 


( 1 ) 


A common method for estimating a calibrator star’s 
size (if it is not known from previous interferometric mea¬ 
surements) is by fitting a photometric energy distrib ution 
(PEP) to m easured photometry (see Appendix [A| . |Boy- 
ajian (2009) found an average difference between angular 


sizes determined by PEP fitting and angular sizes mea¬ 
sured by interferometry to be ~10%, so a 10% error in the 
angular size is adopted for the calibrator stars observed 
for this work. Small calibrator stars are used because 
the smaller a star is, the less its estimated intrinsic visi¬ 
bility is affected by inaccuracies in its size estimate. For 
example, a small calibrator with a 10% error (angular di¬ 
ameter, 0 — 0.2 ± 0.02 mas) observed with the CHARA 
Array’s longest baseline (B = 331 m) in the AT-band will 
have an estimated intrinsic visibility of 0.974 ± 0.005 (a 
0.5% error due to the inaccuracy of an PEP-determined 
size). A calibrator that is twice as large (0 = 0.4 ± 0.04 
mas) and observed in the same way will have an esti¬ 
mated intrinsic visibility of 0.90 ±0.02 (a 2.2% error due 
to the inaccuracy of an PEP-determined size). As a rule 


of thumb, good calibrators are ones that are smaller than 
approximately half the resolution of the observation to 


avoid significant errors in the calibrator’s visibility (van 
Belle fc van Belle||2005 ). 


For this work, at least two calibrator stars were ob¬ 
served for each target star to help mitigate calibrator 
size errors. Their angular diam eters are estimated by 
fitting PHOENIX model PEPs (Husser et al. 2013) to 
photometry gathered from the literature. Three of the 
calibrators used here (HP 177003, HP 185872, and HP 
201908) had temperatures greater than the PHOENIX 
model grid (which goes as high as 12000 K). For these 
three calibrat ors, the PEP fits were m ade using ATLAS9 
model PEPs (Castelli & Kurucz|^04|). Calibrator angu¬ 
lar diameters are listed in 'Tabled and range from 0.156 
to 0.582 mas. 


4. FUNDAMENTAL STELLAR PROPERTIES 
4.1. Oblate Star Model 

The limb-darkened disk model traditionally used to 
analyze interferometric visibilities takes neither the dis¬ 
tended shape of rapidly rotating stars nor the gravity 
darkening caused by this distended shape into account. 
The model used here employs a Roche geometry and is 
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Figure 2. The phot osph eres of the best fitting models for the three UMa nucleus stars modeled here: |2a| - Phecda (HD 103287), |2b| - 
Megrez (HD 106591), |2c| - Alcor (HD 116842). The black points represent a grid of colatitudes and longitudes on the near side of the 
model star. The blue circles represent the a uniform disk radius fitted to each individual visibility at the appropriate baseline orientation 
observed. The data are duplicated at 180° orientation. 


based on the mod els u sed in van Belle J 20121 , Aufden- 
berg et al.| (|2006|), and |Monnier et ai.| (|2U12p. In order 
to determine the fundamental properties ot rapid rota¬ 
tors, the observed visibilities (Section]^ and broadband 
photometry are compared to model-predicted visibilities 
and photometry; the adopted photometry for each star 
is assembled in Appendix]^ The eight input parameters 
for the model star are its equatorial radius (Re), its mass 
(M*), its equatorial rotational velocity (K)) the inclina¬ 
tion of its polar axis relative to our line-of-sight (j), the 
gravity darkening coefficient used in the model (/3), the 
temperature at its pole (Tp), the parallax of the observed 
star (TTpix), and the position angle of its pole ('0) with a 
180° ambiguity. Of these, the parallax is set by Hip- 
parcos measurements, the gravity darkening coefficient 
is set by one of two possible relations (see below), and 
the mass is estimated from evolution models (see below). 
The remaining five parameters (i?o, 14, i, Tp, and ip) are 
allowed to vary under the constraint that the equatorial 
velocity (14) must yield a model vsini that is consistent 
with the observed wsinz. 


Two gravity darkening laws are incorpo rated here¬ 
with the canonical g ravity darkening law (von Zeipel] 


1924a|b Claret 2000 hereafter, vZ), the stars modeled 


here are hot enough to have fully radiative envelopes 
giving them a gravity darkening coefficient, /3, of 0.25. 
However, a modern gravity darkening law, tested with 
results from interferometric observations of r apidly ro¬ 
tating stars (Espinosa Lara & Rieutord 2011, hereafter 
ELR) shows that pi is dependent on the angular rota¬ 
tion rate, w, and ranges from 0.25 for a non-rotating star 
{oj = 0) to ~0.09 for a star rotating at its breakup veloc¬ 
ity (w = 1). 

The oblateness of a star depends not only on its ro¬ 
tation, but also its mass. After the best fitting free pa¬ 
rameters are determined, the age and mass are calculated 
using evolution models. The mass used in the oblate star 
model is then updated to match the mass determined by 
the evolution model. The oblate star model and evolu¬ 
tionary model are run iteratively until neither the mass 
nor the free parameters change by more than ~0.1% af¬ 
ter a series of consecutive runs, corresponding to an 
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- 0.3 - 0.2 - 0.1 0.0 0.1 0.2 0.3 

East - West (mas) 


(c) 

Figure 3. Same as Figure]^ but for the three UMa stream stars modeled: |^- Chow (HD 141003), |3b| - 16 Lyr (HD 177196), [^- 59 Dra 
(HD 180777). The baseline orientations of 16 Lyr and 59 Dra are undersampled, making it difficult tomeasure their oblateness directly. 


of ^0.002 Rq, a T4 of ~0.2 km s“^, an i of ~0.1°, a Tp 
of ~8 K, and a tp oi ~0.4°. 

To determine the ages and masses of the rapidly rotat¬ 
ing stars in this paper the star’s average radius (Ravg), 
total luminosity (Ttot), and equatorial velocity (Re)) as 
determined by the oblate star model are comp ared to 
the predict i ons o f MESA evolutionary models (Paxton 
et al. 2011, 2013). These three parameters {Ravg, Ttc., 
and Ve) correspond to a star with a unique mass, age 
and angular velocity. The mass used by the oblate star 
model is set equal to the mass determined by this com¬ 
parison in the iterative process described above. For this 
project, MESA evolution track^ were computed for a 
grid of masses and angular velocities (with resolution 
of 0.1 Mq and 10% breakup velocity, resp ectively) at 
a met allicity of Z = 0.016 as measured by [King et al. 
(2003) for the UMa moving group. 

T he stellar model is constructed by calculating the stel¬ 
lar intensity at each point on an oblate spheroidal grid, 
constructed of 51 points along the colatitudinal axis (i?) 


^ See the supplemental material for examples of MESA inlists 
used in this project. 


and 51 points along the longitudinal axis {(p) for a total 
of 2601 points on the star. Then, a radius (R(i?)) and 
surface gravity ((?('d), with radial component, 5r('i^) and 
polar co mponent, q,')('d)) are calculated for each point on 
the grid ( van Bellq|2012 ): 


R(i?) = 3 


R 


wsint? 


,p ^^3[ ^ + arccos(t.sind) j 


g{-d) = \/ 5 r(d )2 + 5,5 (^) 2^ where 
9rW = + RWi^smt}^ (3) 

5j?(i?) = i?(i?)n^ sinilcosi?. 


In this prescription, Rp is the model star’s polar radius: 


“ [ D 
Itn 




1-1 


2GM* 


( 4 ) 


uj is the angular velocity of the star relative to its critical 
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Table 3 
Observing Log. 


Target Name/HD 

Cal HD 

Cal Diameter (mas) 

Combiner 

Baseline 

Bandpass 

^ brackets 

visibilities 

Date 

Phecda 

99913 

0.582 ± 0.058 

Classic 

E2-W2 

K 

2 

2 

4/23/2012 

103287 

99913 

0.582 ± 0.058 

CLIMB 

S2-E2-W2 

K 

2 

6 

6/2/2012 


105525 

0.392 ± 0.039 

CLIMB 

Sl-El-Wl 

K 

2 

6 

5/11/2013 


99913 

0.582 ± 0.058 

CLIMB 

Sl-El-Wl 

K 

3 

9 

5/11/2013 

Megrez 

108954 

0.451 ± 0.045 

CLIMB 

Sl-El-Wl 

H 

4 

12 

4/20/2012 

106591 

108845 

0.481 ± 0.048 

CLIMB 

Sl-El-Wl 

H 

2 

6 

4/21/2012 


108954 

0.451 ± 0.045 

CLIMB 

Sl-El-Wl 

H 

2 

6 

4/21/2012 

Alcor 

119024 

0.306 ± 0.031 

CLIMB 

Sl-El-Wl 

H 

4 

12 

4/20/2012 

116842 

108954 

0.451 ± 0.045 

CLIMB 

Sl-El-Wl 

H 

1 

3 

4/21/2012 


118232 

0.465 ± 0.047 

CLIMB 

Sl-El-Wl 

H 

2 

6 

4/21/2012 

Chow 

140160 

0.293 ± 0.029 

CLIMB 

Sl-El-Wl 

H 

2 

6 

4/21/2012 

141003 

137510 

0.525 ± 0.053 

CLIMB 

Sl-El-Wl 

H 

2 

6 

4/21/2012 

16 Lyr 

177003 

0.156 ± 0.016 

PAVO 

S2-E2 

R 

3 

69 

7/10/2012 

177196 

172883 

0.181 ± 0.018 

PAVO 

S2-E2 

R 

2 

46 

7/10/2012 


177003 

0.156 ± 0.016 

PAVO 

E2-W2 

R 

3 

69 

8/4/2013 


185872 

0.256 ± 0.026 

PAVO 

E2-W2 

R 

3 

69 

8/4/2013 


177003 

0.156 ± 0.016 

PAVO 

E1-W2 

R 

3 

69 

8/5/2013 


185872 

0.256 ± 0.026 

PAVO 

El-Wl 

R 

2 

46 

8/5/2013 

59 Dra 

184102 

0.263 ± 0.026 

PAVO 

S2-E2 

R 

3 

69 

7/10/2012 

180777 

201908 

0.187 ± 0.019 

PAVO 

S2-E2 

R 

3 

69 

7/10/2012 


184102 

0.263 ± 0.026 

PAVO 

E2-W2 

R 

3 

69 

8/4/2013 


201908 

0.187 ± 0.019 

PAVO 

E2-W2 

R 

3 

69 

8/4/2013 


velocity, flcrit: 


spectrum of the star, Fy. 


UJ = 



Wo(l - Wo) 


2 


Wo 


Ve^Rp 

2GM*’ 


(5) 


and n is the angular velocity of the star in radians per 
second: 

— Luflcrit — 

This allows the gravity dependent surface temperature 
(r(i?)) to be calculated at each point on the grid: 


pTT p‘2tT 

R>~= / h{'d,>f )Ol^{'d)s\n{'d) fi{'d,(p)d(fd'd (9) 

J-d=0Jip=0 

is the intensity spectrum given by the 
PHOENIX model. is the angular radius of the 

model star as a function of colatitude, (p) is the co¬ 
sine of the angle between the observer and the normal of 
the star: 




—— [—gr{'&)(sm{'d) sin(i) cos(vj) + cos('d) cos(i)) — 

gW 

g^{'d){sin{i) cos{(p) cos(i9) — sin(i?) cos(i))]. 

( 10 ) 


T = r„ 


gW 

gp 


( 7 ) 


where is the surface gravity at the model star’s pole: 


gp = 


GAA 

Ef 


( 8 ) 


A gri (f1of PHOENIX atmosphere models (Husser et al. 
'2013|) are interpolated to determine the intensity spec¬ 
trum at each point on the stellar model surface grid based 
on the temperature and surface gravity of those points. 

Model photometry is calculated by integrating the 2601 
intensity spectra that cover the star to compute the flux 


Note that 95 ) = 0 for /x < 0 (i.e., only light directed 
at the observer is included in the integration). The re¬ 
sulting flux spectrum is convolved with the appropriate 
bandpass filter to compute the specific flux from which 
the photometry is calculated. 

The bolometric flux is simply Eboi = JF\dX and the 
apparent luminosity is then Tapp = 47rTboid^. The to¬ 
tal luminosity, Ttot, is calculated by determining Jx, the 
specific irradiance on each point 

J\{TeS,g)=[ Ix{TeB,g,p)gdg (11) 
integrating over all wavelengths: 


Grid step sizes are 0.5 in logg and 200 K in Tg// 


g) dA 


3 


Tboi(i?) = 27r 


( 12 ) 
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and integrating over the model star’s surface: 


L 


tot 





(13) 


Model visibilities are calculated by first creating an 
image of the model star in the bandpass of the obser¬ 
vations. For example, if the visibilities are observed in 
7J-band, the intensity spectra at the different points in 
the image are convolved with an i/-band hlter. A 2D 
fast Fourier transform (FFT) is taken of that synthetic 
image. This image is 4900x4900 pixels with ~1000 of 
those pixels (in the center of the image) being made up 
of synthetic starlight. This distribution is designed to 
produce an image that is high enough resolution to de¬ 
tect the oblateness and for the FFT to extract accurate 
visibilities. The model squared visibility is the complex 
square of that transform at the observed u and v spatial 
frequencies and the model visibilities are the square root 
of that quantity. 

The above prescription yields visibilities and photome¬ 
try based on a model star that can be tuned to match the 
observations. A random search algorithm is employed 
to find the set of free-parameters (i?o, K) h Tp, and 
tjj) that minimizes the difference between observed and 
model predictions. For each set of input parameters, a 
reduced goodness-of-fit metric is calculated with five 
degrees of freedom for both the visibilities and the pho¬ 
tometry. The final (hereafter, Xtot) then calculated 
by adding the values of the visibility data and those of 
the photometry, assuming equal weight for the two. The 
search algorithm randomly selects a set of parameters 
within a given window of parameter space. The initial 
window size for the set of parameters, (i?c, 14, i, Tp, ip) 
is (±0.5 Ro, km s-4 ±20°, ±500 K, ±30°), and 

this search area is decreased over multiple steps, eventu¬ 
ally reaching (±0.01 Rq, ±1 km s“^, ±1°, ±1 K, ±1°). 
This window is initially centered on the initial guess pa¬ 
rameters, but it is re-centered whenever a model with 
a smaller Xtot is calculated. The best fitting model is 
determined by minimizing the Xtot after multiple itera¬ 
tions. The error for each of the five free-parameters is 
found by first varying them independently until the Xtot 
increases by 1 after first scaling the Xtot such that the 
minimum Xtot = 1- 

Due to the large scatter in the broad-band photomet¬ 
ric measurements relative to their error, the best fitting 
model finds an unsealed Xtot of > 100 (dominated by 
the photometric x^) when adopting the published er¬ 
rors for the photometry measurements, the mean and 
median of which are 0.016 and 0.011 mag, respectively. 
More importantly, few of the photometric measurements 
overlapped with the model RED which could indicate 
underestimates of the photometric error, inaccuracies of 
the synthetic spectral energy distribution, incorrect filter 
profiles or zero-points, etc. To account for this, photo¬ 
metric errors of 0.03 mag were adopted for all photomet¬ 
ric values which had an error less than 0.03 mag. With 
these adopted photometric errors, all of the best fitting 
models had an unsealed Xtot of < 15- 

To determine the errors in the age and mass, the age 
and mass are calculated for the ten points which repre¬ 
sent the Icr-errors of the five parameters in the oblate star 


model (i.e., [Re ± 14, 4 Fp, '4], [Re, K ± cry^, 4 Tp, 

Ip], etc.). The lowest and highest values that come from 
this procedure represent the lower and upper bounds of 
the statistical errors presented here. We note that this 
method does not take into account any correlations that 
may be present between the free parameters. The final 
best fitting parameters and their errors for the vZ and 
ELR gravity darkening laws can be found in Tables!^ and 
[§ respectively. Figures ^ - [9] illustrate the best fitting 
model visibilities and photometry, as constrained by the 
observations, for both gravity darkening prescriptions. 


4.2. Initial Model Parameters 

The x^ minimization technique that is u sed to deter¬ 
mine the best-htting model (see Section 4.1) is especially 
sensitive to the initial guess given for tire star’s inclina¬ 
tion. To account for this, for each star, the model is run 
a number of times using various fixed inclinations. The 
inclinations chosen range from 90° (edge-on) down to an 
inclination that would have the model star rotating at 
breakup velocity given its vsini. The best-htting set of 
parameters of these hxed-inclination models is chosen as 
the set of in put parameters for the process described in 
Section 14.11 

The initial guess value for M* that is supplied for the 
model runs at hxed inclinations is determined based on 
the star’s spe ctral type a nd the spectral type-mass rela¬ 
tions found in Cox (2000). The initial guess values for Ro 
and Tp are based on the angular diameters and effective 
temperatur es listed in the JMM C Stellar Diameter Cata¬ 
log (JSDC, Lafrasse et al.|20l0 ) for each star. The initial 
value for ip is determined by fitting a uniform ellipse to 
the visibilities in the cases where multiple baseline orien¬ 
tations have been used or is set to 0° in the cases where 
they have not. 


4.3. Merak 

The apparent slow rotator, Merak (HD 95418), was ob¬ 
served using the Clas sic beam combiner on t he CHARA 
Array previously by ||Boyajian et al. (2012). We have 
taken the radius and luminosity determined by that 
study as well as its v sin i to determine its age and mass 
using the MESA evolution model using a similar process 
described in Section |4.1[ but without any iteration. Be¬ 
cause of this, we do not determine the inclination of this 
star nor its equatorial velocity. We assume an edge-on 
inclination of 90°. The results are compiled in Table 


4.4. 16 Lyr and 59 Dra 

When running the model described above, the two 
stream stars, 16 Lyr and 59 Dra, both yield best htting 
values for Ravg and Ltot that correspond to unphysical 
positions below the zero-age main sequence for their re¬ 
spective best fit values for 14. One way to reconcile this 
discrepancy would be for the stars to have a metallicity 
of Z < 0.013 (~0.1 dex lower than the moving group). 
We are cautious against advocatiM for this interpreta¬ 
tion since, as discussed in Section we have insufficient 
baseline orientations to fully measure the oblateness and 
gravity darkening in these cases. We note that the best 
fitting values for 14 for both 16 Lyr and 59 Dra are suf¬ 
ficiently large that they shift the zero-age main sequence 
above the best htting values for Ravg and Ltot- If these 
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Ve values are too large, this could explain the unphysical 
i?avg and Ltot without changing the metallicity. Figure 
illustrates how the zero-age main sequence is raised by 
rapid rotation. With this in mind, we run the model for 
these two stars constraining the equatorial velocity to be 
within the more modest range of 94 to 202 km s“^ for 
each star. This range corresponds to the dispersion about 
the maximum of the probability distribution of equato¬ 
rial rotationvelocitiesjor late-type A-stars as determined 
by Zorec & Royer (2012). We make this constraint by 


fixing the stars’ inclinations such that i = arcsin(^^^) 

where E[14] is the maximum of the aforementioned prob¬ 
ability distribution. This corresponds to inclinations of 
^57° and ~28° for 16 Lyr and 59 Dra, respectively. 


5. MODEL RESULTS 
5.1. Photospheric Properties 

Using the procedure described in Section [4T] the best 
fitting models for all six of the stars observed show 
values ranging from 3.1—13.4. The model fitting us¬ 
ing the vZ gravity darkening law yield a high inclina¬ 
tion {i > 70°) for one star (Alcor), moderate inclina¬ 
tions (40° > i < 70°) for two stars (Megrez and Chow), 
and a low inclination {i < 40°) for one star (Phecda); 
both 16 Ly r and 59 Dra have fixed inclinations (see 
Section 4.4). These results also show an oblateness, 
p = {Ra — Rp)/IR that ranges from 3% to 54% with 
an average of 26% and temperature differences across 
the photosphere, AT = Tp — Tg that range from 214 K 
to 6414 K with an average of 2965 K. The same analy¬ 
sis done using the ELR gravity darkening law also yields 
a high inclinations for Alcor, moderate inclinations for 
four Megrez and Chow, and a low inclination for Phecda. 
These results show an oblateness range of 3% to 55% with 
an average of 24% and temperature differences across the 
photosphere that range from 192 K to 3769 K with an 
average of 1696 K. The smaller mean temperature gra¬ 
dient seen with the ELR law is because that law yields 
a smaller gravity darkening coefficient, /3, which lessens 
the effect the local surface gravity has on the local tem¬ 
perature. Using the vZ law, /3 is 0.25 for all four observed 
rapid rotators. The ELR law has (3 ranging from 0.138 
to 0.242. 


is partially dependent on the mass because the radius, lu¬ 
minosity, and temperature of more massive stars evolve 
more rapidly, thus allowing for a more precise determi¬ 
nation of the age because fixed uncertainties in these pa¬ 
rameters will correspond to a smaller percent error in the 
age. We caution that these uncertainties are only statis¬ 
tical. Systematic uncertainties (such as those in gravity 
darkening and metallicity) can lead to more substantial 
errors. Only Chow shows a disparity in its age estimates 
between the two gravity darkening laws. Chow’s age is 
determined to be 659l(() Myr when using the vZ law or 
610^14 

Myr when using the ELR law. The final ages and 
masses for the are presented in Table 

5.3. Comparison with Other Evolution Models 

In order to test the accuracy of the MESA evolution 
models and to begin to address some of the systematic 
errors that may be introduced by them, we compare the 
results from one of the stars in our sample across four dif- 



but without accounting for rotation. We use the total lu¬ 
minosity, average radius, and equatorial rotation velocity 
determined for Alcor (HD 1168421^ as our point of com¬ 
parison between the four models. We chose Alcor for this 
comparison because it is the only rapidly rotating nucleus 
member whose rotation speed is less than the maximum 
predicted by the Geneva models, which are restricted to 
values of w of ;< 0.9 for the masses and ages in question. 
The results are listed in Table ID 
The absolute ages agree extremely well between the 
two rotating models, with a percentage difference of only 
0.5% (0.02-cr). The determined stellar masses also show 
good agreement, with a percentage difference of 3.1% 
(1.4-cr). The ages determined by the non-rotating mod¬ 
els also agree with each other extremely well with a per¬ 
centage difference of 0.9% (0.07-cr), but as expected, they 
are systematically older than those determined using the 
models that account for rotation. The masses determined 
by the non-rotating models also show good agreement 
with each other with a percentage difference of 2 . 1 % ( 1 . 0 - 
cr). 


5.2. Masses and Ages 

The masses calculated by the procedures discussed in 
Sectionj^range from 1.447 to 2.509 Mq for all seven stars 
in the sample using the vZ gravity darkenin g law and 
1.443 to 2.509 Mq using the ELR law (Figure 10). The 
mass estimates for the individual stars are consistent be¬ 
tween the two laws within their 1-3% uncertainties with 
the exception of Chow, whose mass is 2 . 333 l'lQ;Qi [5 
using the vZ law or 2.388l'lQ g 2 i using the ELR law. 
The ages calculated by the procedures described above 
range from 401 to 659 Myr for all seven stars in the 
sample using the vZ gravity darkening law and 333 to 
610 Myr using the ELR law. With the exception of the 
star Chow, these age estimates are consistent with being 
coeval using either the vZ and ELR laws, despite their 
larger uncertainties, that range from 2 to 41% and with 
a mean and median uncertainty of 14% and 12%, respec¬ 
tively. It is worth noting that the uncertainty in the age 


5.4. A New Age Estimate for the UMa Moving Group 

The mean age, uncertainty in the mean, and standard 
deviation of the 7 Ursa Major moving group A-stars pre¬ 
sented here are 451, 32, and 86 Myr when using the vZ 
gravity darkening law and 451, 37, and 98 Myr when us¬ 
ing the ELR law. These large standard deviations are 
due in large part to the relatively old age we estimate for 
Chow (6591 Myr for the vZ law or 6 IOI 35 Myr for the 
ELR law). 

The discrepant age for Chow questions its association 
with the moving group. Of the seven stars studied here, 
Chow is one of two stars considered to be a “probable 


member” by King et al. (2003); the other five a r e clas - 


sified as merhbers. As assembled in [King et ah] ( |2003[ ) 
its space motion is consistent with that of nucleus mem¬ 
bers, despite being 23 pc further away (Table [^. Since 

Using the vZ gravity darkening law 
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Table 4 

Comparing Evolution Models. 

Fundamental Parameters for Alcor (HD 116842) 
Average Radius (R©) 1.846 ^q'q57 

Total Luminosity, Ltot (L©) 

Equatorial Velocity (km s“^) 238.6I*]g*^2'^ 

MESA (with rotation) 

Age (Myr) 422l®| 

Mass (Mq) 1.842+g-“3j 

Geneva (with rotation) 

Age (Myr) 4241®^ 

Mass (Mq) 1.899ig;g^^ 

MESA (without rotation) 

Age (Myr) 5751^^ 

Mass (Mq) l-817to,o27 

Padova (without rotation) 

Age (Myr) 580l®^ 

Mass (Mq) 1.855t|5;o29 


we cannot confidently exclude Chow as a member, we 
give statistics both with and without it. If Chow is ex¬ 
cluded, we determine a mean age and standard deviation 
for the 6 remaining stars to be 416 ±11 Myr when using 
the vZ law and 424 ± 79 Myr when using the ELR law. 

A primary goal of this work is to use the ensemble of 
stellar ages to provide a new, independent age estimate 
for the Ursa Major moving group. The distributions 
of individual ages in Figure however, illustrates the 
challenge of doing this robustly as the determined ages 
contain systematic uncertainties (e.g., gravity darken¬ 
ing), a broad range of statistical uncertainties (that can 
bias wei ghted values ) , and possible non-members (e.g., 
Chow). Beers et al. (19901 discuss a variety of statisti¬ 
cally robust techniques for computing the central loca¬ 
tion (“mean”) and scale (“dispersion”) of small samples 
that are potentially contaminated with outliers or that 
have and unknown underlying distribution. Following 
their recommendations, we choose to compute a median 
for the central location of the age and use a technique 
known as t he “gapper” to estimate the dispersion in our 
sample (see Wainer fc Thissen|1976 1 . A median is better 
in this case because it is influenced much less by any in¬ 
dividual point than a mean would be. A median is also 
preferred over a weighted mean for this sample because of 
the broad range of uncertainties that may not account for 
all systematic uncertainties. The gapper method is based 
on the size of the intervals (or “gaps”) in an ordered set 
of measurements with the “gaps” near the median being 
weighted more heavily. The gapper is normalized such 
that it is equivalent to a standard deviation. The median 
age and gapper scale (ug) of the seven A-stars presented 
here are 415 ± 71 Myr when using the vZ law and 408 
±110 when using the ELR law. 

Since the gapper scale is intended to approximate the 
standard deviation for a Gaussian distribution, we use it 
to define an uncertainty in the median as -f=, following 


standard convention. The median, gapper scale, uncer¬ 
tainty in the median, mean, and standard deviation are 
presented in Table for three distinct subsamples of the 


seven stars observed. The first of these subsamples is the 
four nucleus stars (Merak, Phecda, Megrez, and Alcor) 
which are considered bona fide members of the moving 
group, and so are of greater interest in determining the 
age of the group. We find a median age and gapper scale 
of 415 ± 6 Myr and 404 ± 55 Myr for the vZ and ELR 
laws, respectively. The second of these samples is the full 
sample of seven stars with an age of 415 ± 71 Myr (vZ) 
and 408 ±110 Myr (ELR). The final sample is the full 
sample excluding Chow which, due to its estimated old 
age, may be an interloper. Without Chow, we find a vZ 
age of 415 ± 13 Myr and an ELR age of 404 ± 88 Myr. 

As discussed in Section [5T| the model results using the 
two gravity darkening laws show no considerable differ¬ 
ence for in div idual stars. The vZ law, as illustrated in 
Figures pUfTH does yield more consistent age estimates 
(tTg = 13 Myr) among the observed stars (excluding 
Chow) than the ELR law does (ug = 88 Myr). However, 
given that many of the uncertainties in the individual 
measurements are as large or larger than the dispersion 
in the age estimates, we consider that this may be a sta¬ 
tistical anomaly. Because of this, we hesitate to favor 
one law over the other. 

To estimate the age of the moving group, we combine 
the following into one set of age estimates: the age of 
Mer ak d etermined using the method described in Sec¬ 
tion |4^ the ages of Phecda, Megrez, Alcor, 16 Lyr, and 
59 Dra as determined using the vZ law; and the ages of 
those same five stars as determined using the ELR law. 
This combined set of ages allow us to sample what our 
technique can achieve by accounting for the full spread 
in ages we estimate using two gravity darkening laws. 
With this combined set, we find the median age and un¬ 
certainty in the median of the moving group to be 414 ± 
23 Myr. 


5.5. Model Precision in the Age Estimate for Isolated 

A-Stars 

Under the assumption that these stars are the same 
age, the resulting coeval ages provide validation of not 
only the model presented here, but also the MESA evo¬ 
lution model and the physics assumed therein. The dis¬ 
persion of ages can be used to quantify the precision of 
this technique when applied to isolated adolescent-age A- 
stars. Only three stars (Phecda, Megrez, and Alcor) of 
the observed seven are both considered bona fide nucleus 
members of the moving group and were fully modeled by 
the technique presented in Section 4.1 The median and 
gapper scale of their six age estimates (an age estimate 
using the vZ law and one using the ELR law for each 
star) is 415 ± 40 Myr. We use this scale value to de¬ 
termine a precision in our model of ^10% for stars with 
masses ranging from ^1.8 - 2.4 Mq and at a few hundred 
Myr age. Therefore when using this technique on field 
A-stars we expect an overall uncertainty of 10% in the 
age estimates. 


6. SUMMARY 

We present new interferometric observations for six A- 
type stars in the Ursa Major moving group and new 
age and mass estimates for these six plus one that was 
observed previously. All observations were made using 
the Classic, CLIMB, and PAVO beam combiners on the 
CHARA Array. Four of the observed stars are known 
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to be rapidly rotating with tisini > 170 km causing 
them to be measurably oblate. 

To properly account for this rapid rotation, a model 
was constructed with a Roche geometry based on eight 
parameters: R^, M*, 14, i, /3, Tp, TTpix, and ij}. Visibil¬ 
ities and photometry were calculated for each star us¬ 
ing model-generated images and PEDs, and then com¬ 
pared to measured visibilities and photometry. Five of 
the model parameters {Ro, 14 ) 4 Tp, and '0) were allowed 
to vary, with 14 constrained by i and the measured v sin b 
Age and mass estimates were made for each of the stars 
in this sample by comparing their modeled average ra¬ 
dius, luminosity, and equatorial velocity of to those pa¬ 
rameters determined by MESA evolution models. The 
mass determined by the MESA model was then used in 
the Roche model and this process was repeated until the 
models converged. 

Two different gravity darkening laws were studied. 
Neither law was favored by the interferometric and pho¬ 
tometric data, nor was either law favored by the final 
age estimates. The dispersion in the age estimates was 
significantly smaller for the ages estimated using the vZ 
law than the ELR law. However, because this dispersion 
is of the same order of magnitude as the statistical un¬ 
certainties in the ages, we consider that this may be a 
statistical anomaly. The age estimated for Chow makes 


it older than the moving group as a whole and is thus ex¬ 
cluded as a potential interloper in our final age estimate. 
Because neither gravity darkening law was favored, we 
combined the ages estimated with the vZ and ELR laws 
to determine the overall age of the moving group. 

By determining the ages of these coeval stars, we val¬ 
idate this technique for use on individual field stars and 
determine a model uncertainty of approximately 10% for 
stars with masses ranging from ~1.8 - 2.5 M©. Using the 
independent technique described here, we find the age of 
the Ursa Major moving group to be 414 ± 23 Myr. This 
result is consistent with previous age estimates for the 
Ursa Major moving group. 
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APPENDIX 


PHOTOMETRY 


The photometry us ed to construct ph otometric er iergy distributions are prese nted in Table |A1| The optical pho¬ 


tometry is taken from Mermilliod (2006[ UBV) and Hauck & Mermilliod (1997 uvby). Values in these compilations 


are adopted over others because of the large numbe r of observations that they average to compute final values, and 


because all eight stars are included in these surveys. Hauck & Mermilliod (19971 do not report errors for Stromgren y 
band measurements of Phecda and Megrez. For these stars, an error of (J.Ol mag is assu med, which is consistent with 
the photometric uncertainties of stars of similar brightness in their survey. In addition, Hauck &: Mermilliod (19971 


do not report any uncertainties in Stromgren photometry for 16 Lyr or 59 Dra. Near-infrared [JHK) photometry are 
assembled from various sources and is e ither already in the Johnson photometric system or converted to it. 2MASS 
JHK photometry are listed in Table [AT] but are not adopted for m any of the sample stars bec ause they are saturated. 

Merak (HD 95418) - J- and if-band photometry is adopted from Morel & Magnenat (19781, and errors of 0.05 mag 
are assumed. No i7-band photometry is available. 

Phecda (HD 103287) - J-band photometry is adopted from Morel & Magnenat (1978), and errors of 0.05 mag are 
assumed. AT-band photometry is adopted from Ducati (2002). INo H-band photonietry is available. 

Megrez (HD 106591) - J-, H-, and AT-band photometry is adopted from Aumann & Probst (1991), and the adopted 
uncertainties are the average uncertainties of that survey. 

Alcor (HD 116842) - J-, H-, and AT-band photometry is adopted from Kidger & Martin-Luis (2 00^ a fter converting 
the JHK measurements found there to the Johnson system with the method found in Alonso et al. j 1994). The adopted 
uncertainties are the average uncertainties of that survey. 

Chow (HD 141003) - J-band photometry adopted from Selby et al. (1988) after converting the J„ measurement 
found there to the Johnson system with the method foun d there. 'I'he a dopted uncertainty is the reported uncertainty 
in that conversion. A'-band photometry is adopted from Ducati (2002). No J7-band photometry is available. 

16 Lyr (HD 177196) -While 2MASS Al-band photometry is unsaturated for this star, it is not adopted because 
including it causes the model to find a best fit with a in the visibility approximately double what it is without the 
Al-band value. No J- or i7-band photometry is available. 

59 Dra (HD 177196) - 2MASS H- and AT-band photometry is unsaturated for this star, so it is adopted. No J-band 
photometry is available. 
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Table 5 


Model Results using the vZ gravity darkening law. 


Phecda 

Megrez 

Alcor 

Chow 

16 Lyr 

59 Dra 

HD Number 

103287 

106591 

116842 

141003 

177196 

180777 

Equatorial Radius, Re (R©) 

o.^oo_o 148 

o c;i 9+0.075 
z.oiz_o 076 

2 002"''0-003 
z.uuz_o 067 

4 486"’"0-003 
^.^00_o 082 

i.OD ^_0 023 

1 594+0.033 
i.oz^_0 035 

Equatorial Velocity, Ve (km s~^) 


310.4t®;® 

238.6l^°2° 

307.8lg g 

101.61111 

I3l.8l“;° 

Inclination, 2 (°) 

27.ll®;t 

52.01IJ 

90.o1°9°o 

44.8ll;l 

“56.9l°°;l 

O-OQ 0+20.5 

^o-^_25.7 

Polar Temperature, Tp (K) 

11138li° 

100301}2| 

8985lji! 

IOO 91 I!! 

824211! 

723 ll!! 

Polar Position Angle, (°) 



154.9l“;| 

161 “i+lO'® 

iOl.CS 20.2 

82 5+13.0 
o^-9 15.0 

'^•-‘--56.9 

Gravity Darkening, 0 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

Angular Rotation Rate, oj 

0 QQQ+0 001 

U.yyy_0,o03 

u.yu ^_0 008 

n QQ 15 + 0.020 
U.090_o 019 

0 999 + 0-001 
u.yyy_o 001 

0 404+^'*^^^ 
u.4u^_0 065 

0 580+0-062 
U.09U_o 098 

Polar Radius, Rp (R©) 

0 qpq”I”0.064 

Qg4 

1 921+0044 

i.yzi_o 044 

1 "790+0.050 

1 . ^zo_o.o50 

Q oio'T+0.045 

O.UO/ _o .038 

1 ^ 99 + 0.023 

i.ozz_o 022 

1 455+0-030 

i.400_o 032 

Average Radius, iJavg (R©) 

z.oo/_Q 079 

9 1 47+0.053 

1.8461°;°®^ 

8 479+0.053 

c5.4/y_o 045 

1 643+^’®^^ 
i.O+9_o.o23 

1 48^^+0-031 
i.485_o 033 

Average Diameter, 9a,vg (mas) 

u.yoz_Q 029 

0.808l°;°20 

0.6841°;°“ 

n fisn +0 010 
U.boU_o.oo9 

n 4 ns+o oo 6 
U.4U8_0.006 

0 507+^'^^^ 
U.OU^_0.oii 

Equatorial Temperature, Te (K) 

4724+^^^ 

69091^®! 

755611“® 

00915+634 

00^0_mQ 

80281°! 

6887111° 

usini (km s~^) 

171.0t|°j 

244.61“;! 

238.611°;° 

2l7.0l®;° 

85.111!:! 

(59 n + 36.8 

'J^.o_5g g 

Total Luminosity, Ltot (L©) 

42.37t|jl 

23.001!;2! 

13.98l°;l® 

52.8711;!! 

10.451°;°° 

4.8611°;!°! 

Apparent Luminosity, Lapp (L©) 

66.94t5|5 

9Q 00 + 1-26 

Zv5.oo_i 28 

11.841°;°! 

64.00l!lf 

10.421°;°° 

5.1261°;!!° 

Age (Myr) 


414+35 

422l!l 

6591“ 

40ll!l 

436+^36 

‘*'JU_203 

Mass (M©) 

9 q/iq+0.055 

^•34«o 060 

z.uoz o 033 

1 S49+0-027 

1-o42 o 031 

9 OOO+0.015 
Z.OCiO o 015 

1 799+0.013 

1 - <^"^0.013 

1 447+0-014 

-L-^^'0.015 

Visibility 

7.646 

2.719 

4.498 

0.763 

1.083 

1.488 

Photometry 

5.798 

3.214 

4.021 

2.329 

6.313 

5.100 

Total 

13.45 

5.933 

8.519 

3.092 

7.396 

6.588 


As discussed in Section [4.4| the inclinations of 16 Lyr and 59 Dra were fixed. 


Table 6 


Model Results using the ELR gravity darkening law. 


Phecda 

Megrez 

Alcor 

Chow 

16 Lyr 

59 Dra 

HD Number 

103287 

106591 

116842 

141003 

177196 

180777 

Equatorial Radius, Re (R©) 

q qoc+0.204 
o.oo<J_o 257 

...,+0.074 

z.oii _0 068 

2 .OO 1 I!;!!! 

4 195+0-002 
‘±.j.yo_o 084 

1.6511!;!!! 

i.5i8l!;!°3°3 

Equatorial Velocity, 14 (km s~^) 

386.3l!°j® 

318.91“;! 

234.ll!!;! 

282.21!°.^^ 

101.31!!;! 

100.91!!;! 

Inclination, i (°) 

28.51!;! 

50.011° 

86 . 8 l!l ’3 

50.ll!;! 

<150 9+30.1 
oo.y_24 7 

090 9+34.2 

^o-^_25.9 

Polar Temperature, Tp (K) 

105201!!! 

9550l“! 

87621!!! 

95391!!^ 

8270l!! 

71641!! 

Polar Position Angle, (°) 

18.41!!;! 

.Q q+ 44.4 
^'-'•y 42.6 

154.0l“;! 

159.81“;! 

13 0+25.9 

J^9.U 25 6 

161.2l!,°2 

Gravity Darkening, 0 

n 1 qQ+0.008 
U.i35_Q Qj^g 

0 170 +OOO 7 

U.l /U_o 007 

^•^^'-0.004 

O.I 61 I!;!!! 

U.z^z_o 003 

0 241+0-006 

u.z4i_0 006 

Angular Rotation Rate, oj 

0 999 + 0-001 
u.yyy_o,oo 2 

u.y/4_o 014 

0 097 + 0.026 

U.52/'_o 026 

0.9851!;!!! 

u.4Ui_o 066 

0 417+0-116 

i -0.168 

Polar Radius, Rp (R©) 

2 . 1861 !;!!! 

1 QQq+0.046 

i.8y3_o 045 

1 Y9Q+0.046 

1 . /zy_o 048 

9.u/u_o 059 

i.uuy_o 023 

1 477+0.031 

' -0.031 

Average Radius, Ravg (R©) 

2 . 500 I!;!!! 

9 1 94^+0.051 
z.iz ^_0 048 

1 849+0-053 

J--o^y_0.055 

q 479+0.061 

o.^i/ ^_0 056 

1 030+0-023 

i.uou_o 024 

1 497+0-032 

i.^y /_0 032 

Average Diameter, 0avg (mas) 

0 91 2+0 032 

u.yiz_Q 044 

0 . 800 I!;!!! 

0.6851!;!!! 

0.6781!;!!! 

0 405+0-OO6 

U.4UO_o 006 

0 . 510 I!;!!! 

Equatorial Temperature, Te (K) 

675i1!!!5 

72441!!! 

76301!!8 

69671!!! 

80661 !! 

69721!!! 

usini (km s~^) 

184.51!!;! 

244 2+^^’® 

233.71!!;! 

216.6lg! 

84 9+10-3 
o‘*-i^_30.9 

^+^-43.6 

Total Luminosity, Ltot (L©) 

44.571°;°® 

22.041!;°! 

13 07+0.72 
40.0/ _o 74 

58.171!;“ 

10.451!;!! 

^.yuu_o 292 

Apparent Luminosity, Lapp (L©) 

64.741!;®® 

99 qq + 1.43 

ZO.OO_j^ 20 

11.851!;!! 

61.721!;!° 

10.421!;!® 

r -| -1 q+0.219 

O.llO _0 216 

Age (Myr) 

qqq+43 

>JOO_g 3 

400l!! 

4541!° 

6 IOI!! 

370l°° 

580l!!° 

Mass (M©) 

9 A^ 9+0-053 

4.^iZ o 060 

9 048+0-035 
2.U4»o 030 

1 090 + 0.027 
J-.ozo 0 030 

9 qoQ+0.036 

z.ooo 0 . 02 i 

1 79 K + O.OI 3 

i. /ZO o.oi4 

1.4431°;!!! 

Visibility 

6.897 

2.664 

4.481 

1.080 

1.141 

1.542 

Photometry x^ 

6.045 

4.133 

4.235 

3.835 

6.265 

5.060 

Total 

12.94 

6.797 

8.716 

4.915 

7.406 

6.602 


As discussed in Section [4.4| the inclinations of 16 Lyr and 59 Dra were fixed. 
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Table 7 

Fundamental properties of Merak (HD 95418). 


Value Source 


Radius (R©) 

3.0210 ± 0.0383 

Boyajian et al 

1 

20121 

Temperature (K) 

9193 ± 56 

Boyajian et al 

1 

20121 

Luminosity, Ltot (L©) 

58.46 ± 0.47 

Boyajian et al 

1 

20121 

V sin i (km s~^) 

46 ±2.3 


Royer et al. > 

2007 

1 


Inclination, i (°) 90 Assumed 

Age (Myr) 408 di 6 This work 

Mass (Mq) 2.509 di 0.005 This work 


Table 8 

Ages and Masses for Individual Stars. 


Star 

Mass (M©) 

Age (Myr) 

Name 

vZ law 

ELR Law 

vZ law 

ELR Law 

Merak 

2.509 ±0.005 

408 ±6 

Phecda 

n oA o-l-0.055 

^.345_Q QgQ 

9 Ul 9"I"0 053 

^•^J-^_0.060 


000+43 

Megrez 

z.uuz_Q 033 

2.048_o_o30 

414135 

400133 

Alcor 

1 Q/19+0-027 

1.842_o 03 ^ 

1 qoq-I- 0.027 
l.»48_o 030 

4221®^ 

4541“ 

Chow 

ry qoO”I"0.015 

z.ooo_0 015 

ry qqQ-|-0.036 

2.38»_o_o21 

6591“ 

610lJ^ 

16 Lyr 

1 '7oq”I" 0.013 
i./zz_o 013 

1 ■T9ir-|-0.013 
-*-• '^^-0.014 

401131 

3701“ 

59 Dra 

1 447+0-014 
0.oi5 

1 443+0.015 
015 

43611“ 

5801“3 


Table 9 

Age Estimates and Uncertainties (in Myr) for Various Subsets 





vZ law 

ELR law 


Combined 



n 

n* 

Mean ± a 

Median di C7g 

Mean ± a 

Median di C7g 

Mean ± a 

Median ± cjg 

\/n 

Nucleus Members 

4 

7 

415 ± 5 

415 ± 6 

399 ± 43 

404 ± 55 

407 ± 34 

414 ± 35 

17 

All Members 

7 

13 

451 ± 86 

415 ± 71 

451 ± 98 

408 ±110 

454 ± 95 

415 ± 93 

35 

All excluding Chow 

6 

11 

416 ± 11 

415 ± 13 

424 ± 79 

404 ± 88 

421 ± 59 

414 ± 56 

23 


n is the number of stars in each subset and also corresponds to the number of age estimates in the vZ and ELR subsets, n* is 
the number of age estimates in the combined subsets and corresponds to 2n — 1. 
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Figure 4. Top Left - Visibility measurements (red circles) for Phecda (HD 103287) are compared to the best fit model visibilities (blue 
squares) assuming the ELR prescription for gravity darkening. Dashed lines connect individual model and measured values and solid 
lines are the error bars. Top Right - Photometric measurements (red circles) for Phecda (HD 103287) are compared to the best fit model 
photometry (blue squares) assuming the ELR prescription for gravity darkening. The spectral energy distribution from which the PED is 
calculated is plotted in grey for comparison. Bottom Left - Same as Top Left, but for the vZ gravity darkening law. Bottom Right - Same 
as Top Right, but for the vZ gravity darkening law. 
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Figure 5. Same as Figure]^ but for Megrez (HD 106591). 
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Figure 6. Same as Figure]^ but for Alcor (HD 116842). 
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Figure 7. Same as Figure]^ but for Chow (HD 141003). 
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Figure 8. Same as Figure]^ but for 16 Lyr (HD 177196) 
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Figure 9. Same as Figure]^ but for 59 Dra (HD 180777) 
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(a) (b) 



(c) 

Figure 10. Distr ibution of stella r masses v ersus age for 7 stars in the Ursa Major mov ing group as determined using the vZ gravity 
darkening law l |10a| l, ELR law ( |10b[ l, and both l |10c[ l with the model described in Section [4.1[ The circles are slowly rotating stars (Ve < 170 
km s~^) and the diamonds are rapidly rotating (14 > 170 km s“^). The black points are nucleus members and the white points are s tream 
members. The red point shows the mass and age of the nucleus member, Merak, that was previously observed by|Boyajian et al.|(|2012|l 
and is discussed here in Section |4.3| In some cases, the size of the statistical error bar is smaller than the size of the symbol. I'he dark 
vertical li nes r epresent the median m th e ages, the shaded regions represent the gapper scale (the standard deviation equivalent discussed 
in Section |5.4[ |. The dotted lines in |10c| connect the age and mass estimates from the two different laws. 
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(a) (b) 



(c) 

Figure 11. Same as Figure but excluding Chow. 
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Table Al 

Adopted Photometry. 


Bandpass 

Merak 

Phecda 

Megrez Alcor 

Chow 

16 Lyr 

59 Dra 


HD 95418 

HD 103287 

HD 106591 HD 116842 

HD 141003 

HD 177196 

HD 180777 

Adopted Optical Photometry 

Mermilliod 1991 

Johnson U 

2.349 ± 0.014 

2.451 ± 0.010 

3.460 ± 0.009 4.260 ± 0.008 

3.816 ± 0.023 

5.281 ± 0.018 

5.440 ± 0.016 

Johnson B 

2.346 ± 0.011 

2.437 ± 0.006 

3.389 ± 0.009 4.176 ± 0.007 

3.731 ± 0.013 

5.199 ± 0.015 

5.442 ± 0.016 

Johnson V 

2.366 ± 0.009 

2.437 ± 0.005 

3.312 ± 0.007 4.009 ± 0.006 

3.667 ± 0.010 

5.013 ± 0.014 

5.136 ± 0.014 

Hauck 1997 

Stromgren u 

3.741 ± 0.022 

3.860 ± 0.016 

4.849 ± 0.011 5.620 ± 0.015 

5.261 ± 0.029 

6.619 

6.699 

Stromgren v 

2.501 ± 0.012 

2.587 ± 0.011 

3.572 ± 0.010 4.400 ± 0.012 

3.922 ± 0.019 

5.412 

5.696 

Stromgren b 

2.349 ± 0.006 

2.426 ± 0.010 

3.350 ± 0.010 4.110 ± 0.011 

3.715 ± 0.012 

5.106 

5.324 

Stromgren y 

2.355 ± 0.005 

2.420 ± 0.010 

3.312 ± 0.010 4.014 ± 0.011 

3.670 ± 0.010 

5.000 

5.120 

Infrared Photometry from Literature 

Cutri 2003 

2MASS J 

2.269 ± 0.244 

2.381 ± 0.290 

3.316 ± 0.248 3.291 ± 0.226 

3.440 ± 0.290 

4.776 ± 0.282 

4.338 ± 0.222 

2MASS H 

2.359 ± 0.164 

2.487 ± 0.174 

3.306 ± 0.252 3.295 ± 0.228 

3.539 ± 0.276 

4.578 ± 0.036 

4.264 ± 0.144 

2MASS K 

2.285 ± 0.244 

2.429 ± 0.288 

3.104 ± 0.338 3.145 ± 0.244 

3.546 ± 0.318 

4.505 ± 0.016 

4.313 ± 0.018 

Ducati 2002 

Johnson K 


2.33 ± 0.02 

1.76 ± 0.01 

3.42 ± 0.01 



Morel 1978 

Johnson J 

2.350 

2.400 

3.110 




Johnson K 

2.350 

2.370 

3.090 




Kidger 2003 

Johnson J 



3.674 ± 0.004 




Johnson H 



3.623 ± 0.004 




Johnson K 



3.631 ± 0.004 




Neugebauer 1969 

Johnson K 

2.38 ± 0.06 

2.34 ± 0.10 





Aumann 1991 

Johnson J 



3.13 ± 0.02 




Johnson H 



3.10 ±0.02 




Johnson K 



3.10 ±0.02 




Selby 1988 

Johnson J 



3.77 ± 0.03 

3.52 ± 0.03 



Johnson K 



3.63 ± 0.03 

3.43 ± 0.03 



Adopted Infrared Photometry 

Johnson J 

2.35 ± 0.05 

2.40 ± 0.05 

3.13 ± 0.02 3.674 ± 0.004 

3.52 ± 0.03 



Johnson H 



3.10 ± 0.02 3.623 ± 0.004 




Johnson K 

2.38 ± 0.06 

2.33 ± 0.02 

3.10 ± 0.02 3.631 ± 0.004 

3.42 ± 0.01 



2MASS H 






4.264 ± 0.144 

2MASS K 






4.313 ± 0.03 


Note - The photometric errors listed here are those reported by the respective surveys, 
photometric points with reported errors < 0.03 mag (see Section 4.11. 


We adopt errors of 0.03 mag for all 
































